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EDITORS” PREFACE. 


THE volume called Higher Mathematics, the first edition 
of which was published in 1896, contained eleven chapters by 
eleven authors, each chapter being independent of the others, 
but all supposing the reader to have at least a mathematical 
training equivalent to that given in classical and engineering 
colleges. The publication of that volume is now discontinued 
and the chapters are issued in separate form. In these reissues 
it will generally be found that the monographs are enlarged 
by additional articles or appendices which either amplify the 
former presentation or record recent advances. This plan of 
publication has been arranged in order to meet the demand of 
teachers and the convenience of classes, but it 15 also thought 
that it may prove advantageous to readers in special lines of 
mathematical literature. 

It 1s the intention of the publishers and editors to add other 
monographs to the series from time to time, if the call for the 
same seems to warrant it. Among the topics which are under 
consideration are those of elliptic functions, the theory of num- 
bers, the group theory, the calculus of variations, and non- 
Euclidean geometry; possibly also monographs on branches of 
astronomy, mechanics, and mathematical physics may be included. 
It 1s the hope of the editors that this form of publication may 
tend to promote mathematical study and research over a wider 


field than that which the former volume has occupied. 
December, 1905. 
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AUTHOR’S PREFACE. 


THIs brief sketch of the Harmonic Functions and their use 
in Mathematical Physics was written as a chapter of Merriman 
and Woodward’s Higher Mathematics. It was intended to give 
enough in the way of introduction and illustration to serve as 
a useful part of the equipment of the general mathematical 
student, and at the same time to point out to one specially inter- 
ested in the subject the way to carry on his study and reading 
toward a broad and detailed knowledge of its more difficult 
portions. 

Fourier’s Series, Zonal Harmonics, and Bessel’s Functions of 
the order zero are treated at considerable length, with the inten- 
tion of enabling the reader to use them in actual work in physical 
problems, and to this end several valuable numerical tables 
are included in the text. 


CAMBRIDGE, MASS., December, 1905. 
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HARMONIC FUNCTIONS. 


ART. 1. HISTORY AND DESCRIPTION. 


What is known as the Harmonic Analysis owed its origin 
and development to the study of concrete problems in various 
branches of Mathematical Physics, which however all involved 
the treatment of partial differential equations of the same 
general form. 

The use of Trigonometric Series was first suggested by 
Daniel Bernoulli in 1753 in his researches on the musical 
vibrations of stretched elastic strings, although Bessel’s Func- 
tions had been already (1732) employed by him and by Euler 
in dealing with the-vibrations of a heavy string suspended from 
one end; and Zonal and Spherical Harmonics were introduced 
by Legendre and Laplace in 1782 in dealing with the attrac- 
tion of solids of revolution. 


The analysis was greatly advanced by Fourier in 1812-1824 
in his remarkable work on the Conduction of Heat, and im- 
portant additions have been made by Lamé (1839) and bya 
host of modern investigators. | 

The differential equations treated in the problems which 
have just been enumerated are 


Oy _ por (1) 


Of Ox” 
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for the transverse vibrations of a musical string; 


OY _ 1 Oy 34 
or Wax Tx (2) 


for small transverse vibrations of a uniform heavy string sus- 
pended from one end; 


+ 4 [7 
y? | = O, 
To (3) 
which is Laplace’s equation; and 
Ou 0 Du 0 
St oa Tay T Be (4) 


for the conduction of heat in a homogeneous solid. 


Of these Laplace's equation (3), and (4) of which (3) is a 
special case, are by far the most important, and we shall con- 
cern ourselves mainly with them in this chapter. As to their 
interest to engineers and physicists we quote from an article 
in The Electrician of Jan. 26, 1894, by Professor John Perry: 


“There is a well-known partial differential equation, which is 
the same in problems on heat-conduction, motion of fluids, the 
establishment of electrostatic or electromagnetic potential, certain 
motions of viscous fluid, certain kinds of strain and stress, currents 
in a conductor, vibrations of elastic solids, vibrations of flexible 
strings or elastic membranes, and innumerable other phenomena. 
The equation has always to be solved subject to certain boundary 
or limiting conditions, sometimes as to Space and time, sometimes 
as to space alone, and we know that if we obtain any solution of a 
particular problem, then that isthe true and only solution. Further- 
more, if a solution, say, of a heat-conduction problem is obtained 
by any person, that answer is at once applicable to analogous prob- 
lems in all the other departments of physics. Thus, if Lord Kel-. 
vin draws for us the lines of flow in a simple vortex, he has drawn 
for us the lines of magnetic force about a circular current: if 
Lord Rayleigh calculates for us the resistance of the mouth of an 
organ-pipe, he has also determined the end effect of a bar of iron 
which 1s magnetized; when Mr. Oliver Heaviside shows his match- 
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less skill and familiarity with Bessel’s functions ın solving electro- 
magnetic problems, he 15 solving problems in heat-conductivity or 
the strains in prismatic shafts. How difficult it is to express exactly 
the distribution of strain in a twisted square shaft, for example, and 
yet how easy it is to understand thoroughly when one knows the 
perfect-fluid analogy! How easy, again, it is to imagine the electric 
current density everywhere in a conductor when transmitting alter- 
nating currents when we know Mr. Heaviside s viscous-fluid analogy, 
or even the heat-conduction analogy! 

“ Much has been written about the correlation of the physical 
sciences; but whén we observe how a young man who has worked 
almost altogether at heat problems suddenly shows himself ac- 
quainted with the most difficult investigations in other departments 
of physies, we may say that the true correlation of the physical 
sciences lies in the equation of continuity 


y lea tay Td) 


Qu y = Ou i =) 1 
oy i da? 


In the Theory of the Potential Function in the Attraction 
of Gravitation, and in Electrostatics and Electrodynamics,* 
V in Laplace's equation (3) 15 the value of the Potential Func- 
tion, at any external point (x, y, £) due to any distribution of 
matter or of electricity; in the theory of the Conduction of 
Heat in a homogeneous solid + V is the temperature at any 
point in the solid after the stationary temperatures have been 
established, and in the theory of the irrotational flow of an 
incompressible fluid 1 V is the Velocity Potential Function 
and (3) is known as the equation of continuity. 


If we use spherical coórdinates, (3) takes the form 


- ar sin 6 20 sin? O Jp? _! 


rf o(rV) I a(sin 3 I ə 4 = א‎ 


* See Peirce’s Newtonian Potential Function. Boston. 

+ See Fourier’s Analytic Theory of Heat. London and New York, 1878; 
or Riemann's Partielle Differentialgleichungen. Brunswick. 

+ See Lamb's Hydrodynamics. London and New York, 1895. 
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and if we use cylindrical coórdinates, the form 


OV 
az 


av ria rry 
or Tr or r op 


3 (6) 

In the theory of the Conduction of Heat in a homogene- 
ous solid,* u in equation (4) is the temperature of any point 
(x, y, 2) of the solid at any time /, and a” is a constant deter- 
mined by experiment and depending on the conductivity and: 
the thermal capacity of the solid. 


ART. 2. HOMOGENEOUS LINEAR DIFFERENTIAL EQUATIONS. 


The general solution of a differential equation is the equa- 
tion expressing the most general relation between the primi- 
tive variables which 15 consistent with the given differential 
equation and which does not involve differentials or derivatives. 
A general solution will always contain arbitrary (1.e., undeter- 
mined) constants or arbitrary functions. 

A particular solution of a differential equation 15 a relation 
between the primitive variables which is consistent with the: 
given differential equation, but which is less general than the 
general solution, although included in it. 


Theoretically, every particular solution can be obtained 
from the general solution by substituting in the general solu. 
tion particular values for the arbitrary constants or particular 
functions for the arbitrary functions; but in practice it 1s often 
easy to obtain particular solutions directly from the differential 
equation when it would be difficult or impossible to obtain the 
general solution. 


(a) If a problem requiring for its solution the solving of a 
differential equation is determinate, there must always be given 
in addition to the differential equation enough outside condi. 
tions for the determination of all the arbitrary constants or 
arbitrary functions that enter into the general solution of the 
equation; and in dealing with such a problem, if the differen- 
tial equation can be readily solved the natural method of pro- 
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cedure 15 to obtain its general solution, and then to determine 
the constants or functions by the aid of the given conditions. 


It often happens, however, that the general solution of the 
differential equation in question cannot be obtained, and then, 
since the problem, if determinate, will be solved, if by any 
means a solution of the equation can be found which will also 
satisfy the given outside conditions, it is worth while to try to 
get particular solutions and so to combine them as to form a 
result which shall satisfy the given conditions without ceasing 
to satisfy the differential equation. 


(0) A differential equation is linear when it would be of the 
first degree if the dependent variable and all its derivatives 
were regarded as algebraic unknown quantities. If it is linear 
and contains no term which does not involve the dependent 
variable or one of its derivatives, it is said to be linear and 
homogeneous. 

All the differential equations given in Art. I are linear and 
homogeneous. 


(c) If a value of the dependent variable has been found 
which satisfies a given homogeneous, linear, differential equa- 
tion, the product formed by multiplyine this value by any 
constant will also be a value of the dependent variable which 
will satisfy the equation. 


For if all the terms of the given equation are transposed 
to the first member, the substitution of the first-named value 
must reduce that member to zero; substituting the second 
value is equivalent to multiplying each term of the result of 
the first substitution by the same constant factor, which there- 
fore may be taken out as a factor of the whole first member. 
The remaining factor being zero, the product is zero and the 
equation ;, satisfied. 


(d) If several values of the dependent variable have been 
found each of which satisfies the given differential equation, 
their sum will satisfy the equation; for if the sum of the values 
in question is substituted in the equation, each term of the sum 
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will give rise to aset of terms which must be equal to zero, and 
therefore the sum of these sets must be zero. 


(e) It is generally possible to get by some simple device 
particular solutions of such differential equations as those we 
have collected in Art. 1. The object of this chapter is to find 
methods of so combining these particular solutions as to satisfy 
any given conditions which are consistent with the nature of 
the problem in question. 


This often requires us to be able to develop any given func- 
tion of the variables which enter into the expression of these 
conditions in terms of normal forms suited tothe problem with 
which we happen to be dealing, and suggested by the form of 
particular solution that we are able to obtain for the differential. 
equation. 

These normal forms are frequently sines and cosines, but 
they are often much more complicated functions known as. 
Legendre’s Coefficients, or Zonal Harmonics; Laplace's Coef- 
ficients, or Spherical Harmonics; Bessel’s Functions, or Cylin- 
drical Harmonics; Lamé’s Functions, or Ellipsoidal Har- 
monies; etc. 


ART. 3. PROBLEM IN TRIGONOMETRIC SERIES. 


As an illustration let us consider the following problem: 
A large iron plate z centimeters thick is heated throughout 
to a uniform temperature of 100 degrees centigrade; its faces 
are then suddenly cooled to the temperature zero and are kept 
at that temperature for 5 seconds. What will be the tempera- 
ture of a point in the middle of the plate at the end of that 
time? Given a* = ף0.18‎ in C.G.S. units. 

Take the origin of codrdinates in one face of the plate 
and the axis of X perpendicular to that face, and let א‎ be the 
temperature of any point in the plate ¢ seconds after the cool- 
ing begins. 

We shall suppose the flow of heat to be directly across the 
plate so that at any given time all points in any plane parallel 
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to the faces of the plate will have the same temperature. 


Then w depends upon a single space-coordinate 1; - = O and 


~ = O, and (4), Art. I, reduces to 
2 


u “Y 
= = = )1( 
Obviously, u = 100° when ¿=0, (2) 
/ =  Owhenx—=0, )3( 
and u = Owhenr=7; (4) 


and we need to find a solution of (1) which satisfies the con- 
ditions (2), (3), and (4). 


We shall begin by getting a particular solution of (1), and 
we shall use a device which always succeeds when the equa- 
tion is linear and homogeneous and has constant coefficients. 


Assume * u = ef* tY, where ₪ and y are constants; substi- 
tute in (1) and divide through by œ*tY¥ and we get y = œf’. 
and if this condition is satisfied, u = e®*tY is a solution of (1). 

u = /6*1978" is then a solution of (1) no mattér what the 
value of 6. 


We can modify the form of this solution with advantage. 
Let 8 = ui, t then u = 07 “#*e#* is a solution of (1), as is also 
Na ו‎ 


By (2), Art. 2, 
7 E G ואה‎ al 


Uu — € - = eT" cos Ut (5) 
is a solution, as is also 
hen ו‎ | 
i= on ah TA) = ETP sin ur; (6) 


21 
and jy is entirely arbitrary. 


* This assumption must be regarded as purely tentative. It must be tested 
by substituting in the equation, and 15 justified if it leads to a solution. 


+ The letter 2 will be used to represent Y — 1. 
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By giving different values to u we get different particular 
‘solutions of (1); let us try to so combine them as to satisfy our 
conditions while continuing to satisfy equation (1). 


u = 079%9*% sin ux is zero when x = O for all values of m; it 
is zero when x = v if u is a whole number. If, then, we write 
u equal to a sum of terms of the form Ae~*”" sin mx, where 
ın is an integer, we shall have a solution of (1) (see (d), Art. 2) 
which satisfies (3) and (4). 


Let this solution be 
u = Ae" “sin x + 4,607 sin 22 + 1,077" sin 34 -+..., (7) 
A,, A,, A,,... being undetermined constants. 
When ¢ = 0, (7) reduces to 
u = A, sin x + 4, sin 2x + A,sin3x+.... (8) 


If now it is possible to develop unity into a series of the 
form (8) we have only to substitute the coefficients’ of that 
series each multiplied by 100 for 4,, A,, A,... in (7) to have 
a solution satisfying (1) and all the equations of condition (2), 


(3), and. (4). 
We shall prove later (see Art. 6) that 


dd 1 | 
t= Al sin z +7 sin ו‎ 5r +... 


‘for all values of x between O and zw. Hence our solution is 


OO I 
u = Al eo sin x + en sin 34 + me sin 5r +.. | (9) 
To get the answer of the numerical problem we have only 
to compute the value of u when x = 2 and 1 — 5 seconds. As 


there is no object in going beyond tenths of a degree, four- 
place tables will more than suffice, and no term of (9) beyond 


| | | . 7 
the first will affect the result. Since sin 5 = 1, we have to 


‘compute the numerical value of 
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E where a*=0.1858 and ¢= 5. 
log 0' = 9.2672 — 0 log 400 = 2.6021 
log + = 0.6990 colog z = 9.5059 — IO 
log a't = 9.9662 — 0 colog 057 = 9.5982 — IO 
log log e = 0.6178 — 10 bzs = 1.7062 
log log e” = 9.6040 — 10 
log e”? = 0.4018 / = 50°.8. 


If the breadth of the plate had been c centimeters instead 
of x centimeters it is easy to see that we should have needed 
the development of unity in a series of the form 


. MX . 2x . א‎ 
A, sin — +4, sin = + Ay Sin wee 


Prob. 1. An iron slab 50 centimeters thick 1s heated to the tem- 
perature 100 degrees Centigrade throughout. The faces are then sud- 
denly cooled to zero degrees, and are kept at that temperature for 
סז‎ minutes. Find the temperature of a point in the middle of the 
slab, and of a point 10 centimeters from a face at the end of that 
time. Assume that 


ye I. 37 .. DAX 
1 = (sin - +- z Sin "+ sin sme... from x = otox =c. 
Ans. 847.0| 4 


ART. 4. PROBLEM IN ZONAL HARMONICS. 


As asecond example let us consider the following problem: 
Two equal thin hemispherical shells of radius unity placed 
together to form a spherical surface are separated by a thin 
layer of air. A charge of statical electricity 15 placed upon 
one hemisphere and the other hemisphere is connected with 
the ground, the first hemisphere is then found to be at poten- 
tial 1, the other hemisphere being of course at potential zero. 
At what potential is any point in the “ field of force” due to 
the charge? 

We shall use spherical coordinates and shall let Y be the 
potential required. Then Y must satisfy equation (5), Art. 1. 
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But since from the symmetry of the problem V is obviously 
independent of ¢,1f we take the diameter perpendicular to the 


3 


ap 


plane separating the two conductors as our polar axis, 1S 


zero, and our equation reduces to 


. OV 
rar V) I a(sin 6 ag zo (1) 
ar” sin 6 og 


V is given on the surface of our sphere, hence 


V= when r=1, (2) 
where 6) = 1 ifo > 0 > 7,306 AO) =o if = >> 


Equation (2) and the implied conditions that V is zero at 
an infinite distance and is nowhere infinite are our conditions. 


To find particular solutions of (1) we shall use a method 
which is generally effective. Assume* that V = RO where R 
is a function of 7 but not of 0, and Ois a function of @ but 
not of 7. Substitute in (1) and reduce, and we get 


2 
1 ra*(rR) o r d(sin O (3) 
R dr — @Osin 0 20 


Since the first member of (3) does not contain @and the 
second does not contain 7 and the two members are identically 
equal, each must be equal to a constant. Let us call this 
constant, which is wholly undetermined, 2.(+2 + 1); then 


at) 
asin O 
r dark) — => etico Asin eg) = m(m-+ 1); 
R dr 6 sin 6 ag 
79 R 
whence r- 5 ) mim + 1( =0, (4) 
{. dO 
I d(sin 3 
and 2 o 277 שק ד‎ + 1)0 = 0. )5( 


* See the first foot-note on page 75. qa 
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Equation (4) can be expanded into 
OR aR 
היק‎ + 270 — m(m t+ 1(20 = 0 
and can be solved by elementary methods. Its complete 
solution 5 | 
R = 447% 4+ Brom} (6) 
Equation (5) can be simplified by changing the independ- 
ent variable to x where x = cos 0. It becomes 


d y , -i 

rane — r7 + m(n + 1)0 = 0 (7) 
an equation which has been much studied and which is known 
as Legendre's Equation. 


We shall restrict 77, which is wholly undetermined, to posi- 
tive whole values, and we can-then get particular solutions of 
(7) by the following device: 

Assume* that O can be expressed as a sum or a series of 
terms involving whole powers of x multiplied by constant 
coefficients. 

Let O = 250" and substitute in (7). We get 


Èin(n — Ija, — n(n + 1)a,1* + m(m + 1)a,x"] = 0, (8) 
where the symbol 2 indicates that we are to form all the 
terms we can by taking successive whole numbers for z. 


“Since (8) must be true no matter what the value of +, the 
coefficient of any given power of x, as for instance 1“, must 
vanish. Hence 


(E + 2E F (6 — HE + 1)a, + mln + 1)a, = o, 


mm+1)— kk+1) 
oT oo A 
FEF) % 9) 
If now any set of coefficients satisfying the relation (9) be 
taken, O = 3a,x* will be a solution of (7). 


and Ary, = — 


If Æ = m, then 000, = 0, ap, =0, etc. 


* See the first foot-note on page 475, 32 
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Since it will answer our purpose if we pick out the simplest 
set of coefficients that will obey the condition (9), we can take 
a set including 2. 

Let us rewrite )0( in the form 


QEYB Niya 
677 (m — 000 + k — 1) (10) 

We get from (10), beginning with א‎ = m — 2, 
z = n(n — 1) 

"-%77 2. (2m -- 1)” 

mim — Din — 20m — 3) 

67. = > Sr, ae 

En 2.4.(2m — 12m -- 3) ” 
rest AI E a 


2.4.6.)27/ — 1212 — 31212 — 5) 


If אל‎ is even we see that the set will end with a,; if mm 
is odd, with a.. | 
- 10)74 - 1 
A 
2 .(2712— 1) 
mim — 1()/ — 211 — 
A | 
2.4.(2m — 1)(2m — 3) 
where a,, is entirely arbitrary, is, then, a solution of )7 It is 
found convenient to take a,, equal to 
(212 — 1)\(2m — 3)...1- 
71 ! | 
and it will be shown later that with this value of 00, O= 1 
when + = I. 


O is a function of ע‎ and contains no higher powers of x 
than z”. It is usual to write it as P,,(7). 


We proceed to write out a few values of /,(x) from the 
formula 


(2m — 1()2% — 3)... 1 E = mm — 1) ב‎ 
Emb) = m! úl 2. (211 — I) 


mim — 1()8 — 2()7% — 3) הי‎ ,/ 


+ 2.4.(2m — 1)(2m — 3) 


(11) 
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We have: 
P(x) = 1 or P(cos 0) = 1 | 
Pix) = א‎ or P cos 0) == cos 0, | 


Piz) = 432° — 1( or Picos 0) = 3(3 cos*0 - 1), 
P(x) = 4(52° — 3x) or P,(cos 0) = 4(5 cos"0 — 3 cos 0), 
P(x) = $(352* — 304° + 3) or 
P,(cos 0( = 4(35 cos'0 — 30 cos’G + 3), 
P(r) = 4(632° — 7014” + 154) or 
P(cos 0) = 4(63 cos’ 6 — 70 cosè’ + 15 cos 0). | 


] (12) 


We have obtained O = £,, (+) asa particular solution of (7), 
and O=P,, ¡cos 0) as a particular solution of (5). (+) or 
P, cos 0) is a new function, known as a Legendre’s Coefficient, 
or as a Surface Zonal Harmonic, and occurs as a normal form 
in many important problems. 


V = r"P, cos 0) is a particular solution of (1), and 7”P,,(cos 0) 
is sometimes called a Solid Zonal Harmonic. 


V = A,P (cos 0) + A,rP,(cos 6) + A,rP,(cos 6) 

+ A r’ Pacos 0) +... (13) 
satisfies (1), is not infinite at any point within the sphere, and 
reduces to 
V = A,P,(cos 0) + A,P (cos 0) + A,P,(cos 0) 

+ A,P,(cos O) +... (14) 


when 7 = I. 


age 0)  A,Pcos 0) 


+ (15) 


satisfies (1), is not infinite at any point without the sphere, is 
equal to zero when 7 = «o, and reduces to (14) when 7 = 1. 


Pa Pa + 


A,P,(cos 0) 4 


y* 


If then we can develop /(0) [see eq. (2)] into a series of the 
form (14), we have only to put the coefficients of this series in 
place of the 4,, 4,, A,,... in (13) to get the value of Y for a 
point within the sphere, and in (15) to get the value of V at a 
point without the sphere. 
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We shall see later (Art. 16, Prob. 22) that if /(@) = 1 for 


o < 6<— and f(A) = o for > > 0 > 7, 
13 2 
HO) = כ 5 (0 ולוו‎ P,(cos 6) 


I 
| ETT 2 A P (cos 6) = .. (16) 


Hence our required solution is 


ES fy 
V=>3+¿>Px4cos 6) 3 "> r`P cos 0) 


E eee AE, (cos#)—... (17) 


12 2.4 
at an internal point; and 


P(cos6)— 5 ZA + P(cos 0(‏ 7 יי 


1.31 
rene P(cos@)—... (18) 


et an external point. 


If ץ‎ = A and 8 = O, (17) reduces to 


I 
45 


an y oe 4 Emb = I. 
, since P,,( 1) 


To two decimal places Y = 0.68, and the point» = -, 0 0 
is at potential 0.68. 
If r = 5 and 0 = > (18) and Table I, at the end of this 


chapter, give 


3 I 7 1.3 I 

V=; q 0-7071 z’ — -7-0.1768 +... 5% 
44 

and the point ? = 5, 0 = e is at potential 0.12. 


If the radius of the conductor is a instead of unity; we 


have only to replace r by : in (17) and (18). 


PROBLEM IN BESSEL'S FUNCTIONS, 2] 


Prob. 2. One half the surface of a solid sphere 12 inches in di- 
ameter is kept at the temperature zero and the other half at 100 de- 
grees centigrade until there is no longer any change of temperature 
at any point within the sphere. ¡Required the temperature of the 
center; of any point in the diametral plane separating the hot and 
cold hemispheres ; of points 2 inches from the center and in the 
axis of symmetry ; and of points 3 inches from the center in a di- 
ameter inclined at an angle of 45° to the axis of יש‎ | 

Ans. 50° § 507% 730.9 ; 1 1 22.0. 


1 
t 


A e 0 | - a P tasn- 
1 | 
ART. 5. PROBLEM IN BESSEL'S FUNCTIONS. 


As a last example we shall take the following problem: 
The base and convex surface of a cylinder 2 feet in diameter 
and 2 feet high are kept at the temperature zero, and the upper 
base at 100 degrees centigrade. Find the temperature of a 
point in the axis one foot from the base, and of a point 6 inches 
from the axis and one foot from the base, after the permanent 
state of temperatures has been set up. 


If we use cylindrical codrdinates and take the origin in the 
base we shall have to solve equation (6), Art. I; or, represent- 
ing the temperature by א‎ and observing that from the sym- 
metry of the problem z is independent of ₪ 


2-7 1 88 | Ou | 


EE = 0 I 

ar” + ror 32 (1) 
subject to the conditions 

u=0 when 2=0, (2) 

u = O 2 >> (3) 

u = IOO “ oa, (4) 


Assume u = RZ where X is a function of r only and 2 of 
z only; substitute in (1) and reduce. 


17R,1dR עלקו‎ 


Weset - 77720 . 


The first member of (5) does not contain z; therefore the 
second member cannot. The second member of (5) does not 
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contain 7; therefore the first member cannot. Hence each 
member of (5) is a constant, and we can write (5) 


IPR, א‎ 102 , 
Rd’ rRdr Zd As (0) 


when yp’ is entirely undetermined. 


rig 
ל‎ : == , 
Hence at H VAS (7) 
PR |, idR 
d ika 2 "R=O, | 
- dr’ To e ° (8) 


Equation (7) is easily solved, and its general solution is 


2 = Ae t Bere, or the equivalent form 
2 = 6 cosh (uz) + D sinh (pz). (9) 


We can reduce (8) slightly by letting ur = x, and it becomes 


Ai 1-6 (10) 


ax x AX 


Assume, as in Art. 4, that A can be expressed in terms of 
whole powers of x Let R = 2009" and substitute in (10). 
We get 

[n(n — 1)a,1*7* + na,x*-* + ["ע,0‎ = 0, 


an equation which must be true, no matter what the value of ₪ 


The coefficient of any given power of x, as #*~’, must, then, 
vanish, and 

k(k -- 1)a, + ka, + a,_,= O, 
or ka, + a-,=— 0, 


whence we obtain 0%- = — Fa, (11) 


as the only relation that need be satisfied by the coefficients in 
order that R = 2a,x* shall be a solution of (10). 


If k =0, 4p_,=0, Q&_,= 0, etc. 


We can, then, begin with א‎ = 0 as the lowest subscript. 
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From (11) 008 = -- 


Then a, a = — a, — y 3) Qe — — "3.3 7-9? etc. 


x y x 
— R=0,|1 ק'ע-‎ o to 
where æ, may be taken at pleasure, is a solution of (10), pro- 
vided the series is convergent. 


Take a, = 1, and then א‎ = /,(x) where 


0 8 


; x" x" x x% 
J) = Ita הת‎ )12( 
is a solution of (10). 

T(x) is easily shown to be convergent for all values real or 
imaginary of x, it is a new and important form, and is called a 
Bessel's Function of the zero order, or a Cylindrical Har- 
monic. 

Equation (10) was obtained from (8) by the substitution of 
x = pur; therefore 


0 0 (ury (urf (ur) 
k=J/ (ur) = 1 - בי"‎ PO 


is a solution of (8), no matter what the value of mu; and 
u = J (ur) sinh (uz) and u = / (ur) cosh (uz) are solutions of 
(1). u = J (ur) sinh (uz) satisfies condition (2) whatever the 
value of u. In order that it should satisfy condition (3) ₪ 


must be so taken that ורוה % יי‎ 
JM = 0: J, ₪ = È (13) 


that is, 4 must be a root of the transcendental equation (13). 


It was shown by Fourier that / (4) = 0 has an infinite num- 
ber of real positive roots, any one of which can be obtained to 
any required degree of approximation without serious difh- 
culty. Let u,, M, M,, - - be these roots; then 


u = A, J hmr) sinh (4,2) + A, J (u,7) sinh (1,2) 
+ A, J (57) sinh (14,2) ape הררה‎ (14) 


is a solution of (1) which satisfies (2) and (3). 
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If now we can develop unity into a series of the form 


t= BS (ur) + B, Jdr) + BL LAI) + 


B, sinh (4,2) ye, pp Ba sinh (ae) roy yt . | (5) 


“= | sinh -sinh (24) sinh (21,) 


satisfies (1) and the conditions (2), (3), and (4). 


We shall see later (Art. 21) that if / (x) = — E 
X 


„| Zr) y JAM?) y Sr) + l 
Hi J, J Aa) ANA H, J (H) 0 


for values of ? > I. 
Hence 


(16) 


J (4,7) sinh (4,2) J ur) sinh (4,2) 
p, J (p) sinh (24) * 4, /,(44) sinh (244) ג‎ 1 (17) 


is Our required solution. 


283 = 200 


At the point 7 = O, 2 = 1 (17) reduces to 


Y = 200| > sinh 0 > Sinh Ha -+ 
0 H, J (1) sinh (24,) | M4, (14, sinh (24) "= | 


I I 
= 100) ->_—— + —____-_____ 4+... |, 
Fy, (4) cosh y, UY (a) cosh pa, Y \ 
since /,(0) = I and sinh (2%) = 2 sinh x cosh x. 


If we use a table of Hyperbolic functions* and Tables II 
and III, at the end of this chapter, the computation of the 
value of u is easy. We: have 


H, = 2.405 = M= 5.520 
J A) = 0.5190 J (a) = — 0.3402 
colog pm, = 0.018( — 10 colog w= 9.2581 — IO 
Ñ VAVA - 0.2848 / Ju.) — 0.408 37 


“cosh 1, = 9.2530 — 10  “coshu= 7.9037 — 0 


9.1567 — IO 7.030In — 10 


* See Chapter IV, pp. 162, 163, for a four-place table on hyperbolic func- 
tions. 


PROBLEM IN BESSEL'S FUNCTIONS, 25 


(1, 7,(u,) cosh (7 = 0.1434 
(1, J (14) cosh (7 = — 0.0058 


0.1370; u = 08 
At the point r = 4, z = 1, (17), reduces to 


0 TACTA, J (Fy) | 
₪ 100| 0) coa u, J (4) cosh y, * pa J (Ha) coshp, > 1 
./,)30,( = 0.6698 


log /,(44,) = 9.8259 — 10 
colog pu, / (4) cosh u, = 0.1507 — IO 


8.9826 — IO; 
J3u,) = — 0.1678 
log /.(34,) = 9.22482 — 0 
colog u,/,(u,) cosh y, = 7.63017% — 0 
6.8549 — 10; 
1 
MJ aE מו‎ 0.0901 
JEM) __ 0.0007 


447, (1,) cosh 4, 0.0968”  u=9"7 


If the radius of the cylinder is a and the altitude J, we have 
' only to replace u by pa in (13); 24,, 2u,,... in the denomi- 
nators of (15) and (17) by 4,0, 6, - ..; and Mis Ma, May + +. IM 
the denominators of (16) and (17) by ma, p,a, u,a, 


Prob. 3. One base and the convex surface of a cylinder 20 cen- 
timeters in diameter and 3o centimeters high are kept at zero tem- 
perature and the other base at 100 degrees Centigrade. Find the 
temperature of a point in the axis and 20 centimeters from the cold 
base, and of a point 5 centimeters from the axis and 20 centimeters 
from the cold base after the temperatures have ceased to change. 


Ans. 13°.9; 6 
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ART. 6. THE SINE SERIES. 


As we have seen in Art. 3, it is sometimes important to be 
able to express a given function of a variable, x, in terms of sines 
of multiples of +. The problem in its general form was first 
solved by Fourier in his “Théorie Analytique de la Chaleur” 
(1822), and its solution plays an important part in most branches 
of Mathematical Physics. 


Let us endeavor to so develop a given function of x, Ax), 
in terms of sin y, sin 27, sin 3x, etc., that the function and the 
series shall be equal for all values of x between O and z. 


We can of course determine the coefficients 4,, @,, @,,... Ga 
so that the equation 


fx) =a sin xy + a, sin 2x + a, sin 39 +... + a, sin nx (1) 
shall hold good for any z arbitrarily chosen values of x between 
O and z; for we have only to substitute those values in turn 
in (1) to get z equations of the first degree, in which the z co- 
efficients are the only unknown quantities. 


For instance, we can take the 7 equidistant values 41, 24x, 


34x,...nÁx, where Ax = , and substitute them for x in 


(1). We get 
/)22( = a sin 49 -ך‎ a,sin24dr+ a,sin34r-+... ) 
+ a, sin zx, 
/0242( =a, sin 24r + a,sin44r+ a,sin64r4+... 
+ a, sin 224, | 
N34r)=a, sin 342 + a, sin 64% - a, sin ọ4r +... (2) 
/ + a, sin 314, 


7 
nt 1 


Kunda) = a, sin 24x +a, sin 224x +a, sin 324x +... 
+ a, sin 2°42, 
n equations of the first degree, to determine the 2 coefficients 
60700 Cay Aggy oo o Aye 
Not only can equations (2) be solved in theory, but they 
can be actually solved in any given case by a very simple and 
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ingenious method due to Lagrange,* and any coefficient a,, can 
be expressed in the form 


Ly = 


2 n> J(K4x) sin (kmx). (3) 

If now א‎ is indefinitely increased the values of x for which 
(1) holds good will come nearer and nearer to forming a con- 
tinuous set; and the limiting value approached by a, will 
probably be the corresponding coefficient in the series required 
to represent /)2( for all values of x between zero and zx. 

Remembering that (z + 1)4x = xx, the limiting value in 
question is easily seen to be 


an = 2% / Ax) sin mrda. (4) 


This value can be obtained from equations (2) by the fol- 
lowing device without first solving the equations: 


Let us multiply each equation -in (2) by the product of 41 
and the coefficient of a,, in the equation in question, add the 
equations, and find the limiting form of the resulting equation 
as m increases indefinitely. 


The coefficient of any a, a, in the resulting equation 8 

sin kx sin 7249. Ax -+ sin 2k4x sin 2 . 49 +.. 
+ sin 412 sin 77/42 . Ax. 

Its limiting value, since (z+ 1(42 = is 


T 


J sin kx sin mz. dax: 
0 
but 


/ sin 2 sin mx . dx = 3 [ [cos (m — K)x — coslin + r)r]dr =0 
0 0 


if 2 and א‎ are not equal. 


* See Riemann's Partielle Differcntialgleichungen, or Byerly’s Fourier’s 
Series ánd Spherical Harmonics. 
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The coefficient of a,, is 
4x(sin" mdx + sin? 2/4/29 + sin? 3mMmdr +... 4 sin? nm4x). 


Its limiting value is 


TT 
ו‎ TE 
/ sin? mx. dx = 5 


0 


The first member is 


(42) sin mdx. Ax + (24x) sin 2/4 . 4 +... 
+ (1.4x) sin ?/%4 . Ax, 


and its limiting value is 


/ fix) sin mx .dx. 


Hence the limiting form approached by the final equation 


as 2 1S increased 1s 


ד 
2 
Ke) sin mx .dr = 5 om‏ / 
0 


Whence i = - fx) sin mx . dx (5) 
\ 0 


as before. 
This method is practically the same as multiplying the 
equation 
Ke) = 0, sin x+ a, sin 24 + a, sin 3x-+... (6) 
by sin mx. dx and integrating both members from zero to א‎ 


It is important to realize that the considerations given in 
this article are in no sense a demonstration, but merely estab- 
lish a probability. 

An elaborate investigation * into the validity of the develop- 
ment, for which we have not space, entirely confirms the results 
formulated above, provided that between x = 0 and x == the 


* See Art. to for a discussion of this question. 
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function is finite and single-valued, and has not an infinite num- 
ber of discontinuities or of maxima or minima. 


It is to be noted that the curve represented by y= f(x) 
need not follow the same mathematical law throughout its 
length, but may be made up of portions of entirely different 
curves. For example, a broken line or a locus consisting of 
finite parts of several different and disconnected straight lines 
can be represented perfectly well by y = a sine series. 


As an example of the application of formula (5) let us take 
the development of unity. 


Here Ni E 
TT 
2 
= — ¿sinmx.dx; 
1A 
0 
>. , COS mx 
sin mx. dx = — ———. 
it 


2 1 1 = 
fain mx.ax = 21 — COS mn) = zL — )- 1) | 


0 


= 0 If 2 is even ` 


2 4 : 
= — if m is odd. 
2 


sin x , sin 3% , sinSx , sin 7x 
e SSE SE g 


Hence I = 1 
TEN I 3 5 7 


50 


It is to be noticed that (7) gives at once a sine development 
for any constant c. It is, 


_4c(sinz | sin3r | sin5x 
ב‎ - T + l +...) (8) 


Prob. 4. Show that for values of x between zero and 7 


sin x sin 2x , sin2x sin 4x 
הרו ו‎ 
I 2 3 4 


₪ א‎ a=. x -H4 mE _ E e 1 
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if f(x) = ₪ foro > ₪ > 7, and f(x) = ₪ — x for= > s > T. 


(c) f(x) = 


2 E x , 2sin2x | sin3x | sin5x , 2sin6x , sin 7x 


הו ; + 6 i =p‏ ד---7 . + "וה = 


7t TT 
if f(x) =1foro<x< 4 and f(x) = o for — > 2 > 7m, 
(7) sinh x = 
=e | L sing —2 sin ax + sin 3x — = sin gx +... |. 


(e) xX = 


2 (7 4) si יא‎ (Es EN 
a G 3 sin x . sin 2x + F 5 sin 3x dia iil Mca 


ART. 7. THE COSINE SERIES. 


Let us now try to develop a given function of v in a series 
of cosines, using the method suggested by the last article. 
Assume 


f(x) = 6 + b, cos x 4 b, cos 2x + b, cos 3x+>... (1) 


To determine any coefficient 6, multiply (1) by cos mx . dx 
and integrate each term from o to 7. 


n 
fe COS mx . dx =O. 
0 , 
T 
fh cos $x cos mx .dx=0, if mand £ are not equal. 
0 
T 
, 4 ae 
$, cos mx ax = —O,, if wis not zero. 
2 
0 


2 pP 
Hence 6, = — J fx) COS 4. dX,’ )2( 
0 


Af #2 is not zero. 
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To get 6, multiply (1) by dx and integrate from zero to x. 


/ 002 = bn, 


fe cos 22.09 = 
0 


1 T 
Hence 2 = q JN (3) 


which is just half the value that would be given by formula (2) 
if zero were substituted for 7. 


To save a separate formula (1) is usually written 


f(x) = $b, + ù, cos x + 6, cos 2x + 0,605 3x +. (4) 


and then the formula (2) will give 6, as well as the other coef- 
ficients. 


Prob. 5. Show`that for valués of x between o and 7 


_T 4/cosx | 005348 | Cos 5x ); 
(a) א‎ = - (+ 3 24 5” ad 


(3) Ax) = = — 


8 = 2x . cos6x cos 10x 
x 


2 ע.‎ 7-0 t 10° + / 


TI 
if f(x) א > 2100ב‎ > >> and f(x) = x —xfor— <a £7; 
2 


605 * COS 3X | COS5x 
AO EE 
(c) Fx) i : 
if Ax) = 1 foro < æ > 7 and f(x) = o for © > ₪ > m, 

(d) sinh æ =2| ¥(cosh z — 1) — (cosh x +1) cos x 


+ ~ (cosh 7 — 1) cos 2x — — (cosh 7 -1( 60539 + .. 1 


0 7-7 2 2 מ 0 2 


I 2 3 4 


Oom COS% 0605822  COS3% COS 
le) x* = — - ES = p == e+...) 
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ART. 8. FOURIER'S SERIES. 


Since a sine series is an odd function of א‎ the development 
of an odd function of x in such a series must hold good from 
x = — 7 to x = xz, except perhaps for the value x = O, where 
it 15 easily seen that the series 15 necessarily zero, no matter 
what the value of the function. In like manner we see that 
if f(x) is an even function of x its development in a cosine 
series must be valid from xr = — wtor=7, 


Any function of x can be developed into a Trigonometric 
series to which it 1s equal for all values of + between — 7 and z. 

Let f(x) be the given function of x. It can be expressed 
as the sum of an even function of ע‎ and an odd function of x 
by the following device: 


09-00-28 6 


identically; but 


K+) + A— +) 
2 


is not changed by reversing 
the sign of x and is therefore an even function of x; and when 


we reverse the sign of x, e is affected only to the 
extent of having its sign reversed, and is consequently an odd 
function of x. 


Therefore for all values of x between — 7 and zx 


Aa) thm 5 = 10 0 cos x + 4, cos 22 + å, cos 34 + ... 


where cos mx .dx; 


_2 Phx) +A +) 
E ie 


an 


d Ha) — A— x) = 
2 


= q, sin 1 +4, sin 2x +- a, sin 37+... 


where Mi == A pi iitim sin mx . dz. 
T 3 2 
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b,,and a, can be simplified a little. 


/ -2 OE 
4 2 


COS MX. dæ 


ES LL ל‎ COS MX . dry ff) 605 MX. dz |; 


but if we replace x by — x, we get 


0 ד- זז 
K- x)cos mx .dx=— f| Ax)cosmx.dr= | flx)cos mx.dx,‏ 
/ / / 


I 
and we have 00 = — Ax) 605 mx. ax. 
In the same way we can reduce the value of a,, to 


f ,‏ ן 
fo sin mx .dx,‏ — 
Hence‏ 
Na) => 6, +6, cos 2 + 6, cos 2x + 6, cos 3% +...‏ 
+a sin z + a, sin 2x + a, sin 31 +..., (2)‏ 


If 
where 20 = 5 cos mx. dx, (3) 
IP... 
and A, = — f fiz) sin mx. dx, (4) 


and this development holds for all values of + between - x 
and 7. 


The second member of (2) is known as a Fourier's Series. 

The developments of Arts. 9 and 7 are special cases of 
development in Fourier's Series. 

Prob. 6. Show that for all values of x from — x to x 


2 sinh 7 


c= ----- a E — 7005 «+ Ecos za — -i cos 3x4 7 cos 4r.. 1 
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2 ]אק גמן5‎ : . 2. 
— —sinx — sin 2x + sin 3x — 4 sin 4x +... \ 
Tí 2 5 Io 17 


Prob. 7. Show that formula (2), Art. 8, can be written 
Ax) = 2 cosf, + c, cos (x — f,) + ₪ cos (2x — p.) 
+e cos (3x — f,)+..., 


where Cm = )6 +0)? and Jum = tan"! // 


m 


Prob. 8. Show that formula (2), Art. 8, can be written 
f(x) = >a sin 6, +e, sin (x + B,) + 4 sin ) + $) 


+ ¢, sin )3 + f,) +. לי.‎ 
where Cm = (am +0, )* and fn = tan, 


m 


ART. 9. EXTENSION OF FOURIER'S SERIES. 


In developing a function of x into a Trigonometric Series it 
is often inconvenient to be held within the narrow boundaries 
x = ל‎ and z = m. Let us see if we cannot widen them. 

Let it be required to develop a function of x into a 
Trigonometric Series which shall be equal to f(x) for all values 
of x between x = — c and x = cC. 

Introduce a new variable 


7 


יי ₪2 ] 


| 0 
which is equal to — x when ע‎ = — 6 and to ל‎ when 2 = 6 


Ke) = f(s) can be developed in terms of 2 by Art. 8, 
(2), (3), and (4). We have 


F(E) =E, cos 2 + 6, cos 22 + 6, cos 32 +... 


+ a, sin 2 + a, sin 22 + a, sin 38 +..., (1) 


where 0 = — ff = COS WZ. dz, (2) 
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and En (Ea) sin 222 . dz, (3) 


1 
— 


and (1) holds good from 28 = — x to 2= T. 
Replace 4 by its value in terms of x and (1) becomes 


2 TX 
6 


I TA דג‎ 
f(x) = 5 % + 4 cos —— + 4, cos + 4, cos He... 


- +...; (4) 


. AA . 7% À 
+ q, sin — t a, sin ד‎ sin? 


and (2) and (3) can be transformed into 


I 7 mux 
00 = — f A) cos — ax, (5) 
I f , IMAX 
> — f A) sin ----- (6) 
and (4) holds good from ע‎ = — cto xr =C. 


In the formulas just obtained ¢ may have as great a value 
as we please so that we can obtain a Trigonometric Series for 
f(x) that will be equal to the given function through as great 
an interval as we may choose to take. 


It can be shown that if this interval 6 is increased indefi- 
nitely the series will approach as its limiting form the double 


integral — Agar f cos a(A — x)da, which is known as a 
— 00 0 
Fourier's Integral. So that 


i /‏ ןז 
Ae) => f fé f cos (A — aa (7)‏ 


for all values of x. 


For the treatment of Fourier’s Integral and for examples 
of its use in Mathematical Physics the student is referred to 
Riemann’s Partielle Differentialgleichungen, to Schlémilch’s 
Höhere Analysis, and to Byerly"s Fourier's Series and 
Spherical Harmonics. 
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Prob. 9. Show that formula (4), Art. 9, can be written 
(x) = - cos f, + ¢, Cos =- A, | + ¢, cos (ie — p, 
| + e, cos ($ — 8, +..., 
where C = (Gy, bm)? and Pn = tan”? 2* 


m 


Prob. 10. Show that formula (4), Art. 9, can be written 
f(x) = 2 sin 6, + 6, sin ₪ + 2 + c, sin (2 =- + A 


+ [+ 


+ 4 sin 5 


where Em = (a, + Om)? and £, = tan”? 2n, 


mi 


ART. 10. DIRICHLET’S CONDITIONS. 


In determining the coefficients of the Fourier’s Series rep- 
resenting f(x) we have virtually assumed, first, that a series of 
the required form and equal to f(x) exists; and second, that 
it is uniformly convergent; and consequently we must regard 
the results obtained as only provisionally established. 


It is, however, possible to prove rigorously that if f(x) is 
finite and single-valued from א — = ע‎ to x = x and has not 
an infinite number of (finite) discontinuities, or of maxima or 
minima between ע‎ = — x and x = 7, the Fourier’s Series of 
(2), Art. 8, and that Fouriers Series only, is equal to f(x) 
for all values of x between — x and 7, excepting the values of 
x corresponding to the discontinuities of f(x), and the values 
zx and — 7; and that if c is a value of x corresponding toa 
discontinuity of f(x), the value of the series when: = ¢ is 


> imit A +€e)+f(c—e)]; and that when x = x or 


x = — 7 the value of the series is 4[./(7) + f (— 7). 
This proof was first given by Dirichlet in 1829, and may be 


found in readable form in Riemann’s Partielle Differential. 
gleichungen and in Picard’s Traité d'Analyse, Vol. T. 


DIRICHLET'S CONDITIONS. 31 


A good deal of light is thrown on the peculiarities of trigo- 
nometric series by the attempt to construct approximately the 
curves corresponding to them. 


If we construct y = a, sin x and y = 4, sin 2x and add the 
ordinates of the points having the same abscissas, we shall ob- 
tain points on the curve 


y = a, sin x + & Sin 2x. 
If now we construct y = a, sin 3x and add the ordinates to 
those of y = 4, sin x + a, sin 22 we shall get the curve 


y =a sin X + 6 sin 2x + a, sin 34. 
By continuing this process we get successive approximations to 


y =a sin 4 + a, sin 2x + a, sin 3x + a, sin 4t ++." 


Let us apply this method to the series 
y= sin x + 4 sin 9 + 1 sin 5r + (1) (See (7), Art. ©.) 


TT 
y =owhen x =0, 7 fromx=0to + =,and o when z = T, 


It must be borne in mind that our curve is periodic, hav- 
ing the period 2x, and is symmetrical with respect to the 
origin. 


The preceding figures represent the first four approxima- 
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tion to this curve. In each figure the curve y = the series, 
and the approximations in question are drawn in continuous 
lines, and the preceding approximation and the curve corre- 
sponding to the term to be added are drawn in dotted lines, 


Prob. 11. Construct successive approximations to the series 
given in the examples at the end of Art. 6. 
Prob. 12. Construct successive approximations to the Maclaurin's 
3 


Series for sinh x, namely א‎ + + 7 +... 


ART. 11. APPLICATIONS OF TRIGONOMETRIC SERIES. | 


(a) Three edges of a rectangular plate of tinfoil are kept 
at potential zero, and the fourth at potential 1. At what po- 
tential is any point in the plate ? 


Here we have to solve Laplace’s Equation (3), Art. I, 
which, since the problem is two-dimensional, reduces to 


ay SV 
i MA RT I 
dx” + oy’ O, ( ) 
subject to the conditions Y = O when x = O, (2) 
V=0O * S=4 (3) 
V=0 “ y = O, (4) 
V=1 * y=0. (5). 


Working as in Art. 3, we readily get sinh fy sin fx, 
sinh By cos fx, cosh By sin Bx, and cosh By cos fx as particu- 
lar values of V satisfying (1). 


INNY 


V = sinh wa sin == satisfies (1), (2), (3), and (4). 


V = + a a sin — +++. 
| 7 E3 ga a 0 K 
q 07 


is the required solution, for it reduces to 1 when y = >. See 


(7), Art. 6. 


APPLICATIONS OF TRIGONOMETRIC SERIES. 39, 

(6) A harp-string is initially distorted into a given plane 
curve and then released; find its motion. 

The differential equation for the small transverse vibrations. 


of a stretched elastic string 5 


OV _ pI 
at Ox"? 00 


as stated in Art. 1. Our conditions if we take one end of 
the string as origin are 


y = 0 when += 0 (2). 
הח 7 פחק‎ (3) 
OV _ e 

ry, = O / = 0, (4). 
7 = א‎ € 2 = O. (5). 


Using the method of Art. 3, we easily get as particular solutions 
of (1) 

Y = sin fx sin aft, = sin fx cos aft, 

y = cos fx sinaft, and y = cos fx cos + 


na 1 
y = sin = cos 2 satisfies (1), (2), (3), and (4). 


Y = > ay Sin == cos = (6) 
_ 2 . MAX 
where i, = 7 | Fo sin — at (7) 


15. our required solution ; for it reduces to fx) when¢=0. See 
Art. o. 


Prob. 13. Three edges of a square sheet of tinfoil are kept at 
potential zero, and the fourth at potential unity ; at what potential 
is the centre of the sheet? Ans. 0.25. 


Prob. 14. Two opposite edges of a square sheet of tinfoil are 
kept at potential zero, and the other two at potential unity; at 
what potential is the centre of the sheet? Ans. 0.5. 


Prob. 15. Two adjacent edges of a square sheet of tinfoil are 


AO HARMONIC FUNCTIONS, 
kept at potential zero, and the other two at potential unity. At 
what potential is the centre of the sheet? Ans. 0.5. 
Prob. 16. Show that if a point whose distance from the end of a 
. . I a | 
harp-string 1s -th the length of the string is drawn aside by the 


0137675 finger to a distance 4 from its position of equilibrium and 
then released, the form of the vibrating string at any instant is given 
by the equation 


m=00 
2bn' > I. MI . MIX marat 
= — sin — sin-—— cos \ 
, (n — 1(7' 20 n / / 


m=i 


Show from this that all the harmonics of the fundamental note of 
the string which correspond to forms of vibration having nodes at 
the point drawn aside by the finger will be wanting in the complex 
note actually sounded. 


Prob. 17.* An iron slab 10 centimeters thick is placed between and 
in contact with two other iron slabs each 10 centimeters thick. The 
temperature of the middle slab is at first 100 degrees Centigrade 
throughout, and of the outside slabs zero throughout. The outer 
faces of the outside slabs are kept at the temperature zero. Re- 
quired the temperature of a point in the middle of the middle slab 
fifteen minutes after the slabs have been placed in contact. 
Given a? = 0.185 in C.G.S. units. Ans. 10°.3. 


Prob. 18.* Two iron slabs each 20 centimeters thick, one of which 
15 at the temperature zero and the other at roo degrees Centigrade 
throughout, are placed together face to face, and their outer faces 
are kept at the temperature zero. Find the temperature of a point 
in their common face and of points ro centimeters from the com- 
mon face fifteen minutes after the slabs have been put together. 
Ans. 227.8 ; 157.1 ; 24 


ART. 12. PROPERTIES OF ZONAL HARMONICS. 


In Art. 4, 2 = P,,(x) was obtained as a particular solution of 
Legendre's Equation [(7), Art. 4] by the device of assuming 
that z could be expressed as a sum or a series of terms of 
the form a,x” and then determining the coefficients. We 


* See Art. 3. 
+ The student should review Art. 4 before beginning this article. 
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can, however, obtain a particular solution of Legendre’s equa- 
tion by an entirely different method. 
The potential function for any point (4, y, 2) due to a unit 
of mass concentrated at a given point (x,, ,מ‎ 2,) is 
יי‎ - A (1) 
Van (7-0 -2( 6 - 4) 
and this must be a particular solution of Laplace’s Equation 
[(3), Art. 1], as is easily verified by direct substitution. 
If we transform (1) to spherical coordinates we get : 


I 


V = A A SS Sc PA ESSENSE 
Vr? — 2rr,[cos 0 cos 0, + sin O sin 9, cos (@—¢,)] + 7; 


(2) 


as a solution of Laplace’s Equation in Spherical Coordinates 


{(5), Art. 1]. 
If the given point (x,, y,, 2,) is taken on the axis of X, as it 
must be in order that (2) may be independent of ¢, 0, =0, and 


1 


Vr -- 2rr,cos 0 +1 (3) 
is a solution of equation (1), Art. 4. 
Equation (3) can be written 
I 7 71 \ + 
= = — 2— 0 a : 
ץ‎ 4 E cos + (4) 


2 


-$ 
and if ? 15 less than 7, ( — 2” cos 6 + | can be developed 


j 
y Ai 


| i y” l : 
Into a convergent power series. Let 27,— be this series, 
7 


1 
m 


2 being of course a function of*@ Then V = = Dn is a 


mi 
1 1 


solution of (1), Art. 4. 

Substituting this value of V in the equation, and remem- 
bering that the result must be identically true, we get after a 
slight reduction 


I d| . dp, 
min + ont HSI sin jin =0; 
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but, as we have seen, the substitution of x = Cos 0 reduces this 
to Legendre’s equation [(7), Art. 4]. Hence we infer that the 
coefficient of the mth power of 2 in the development of 
)1 — 222 -ך‎ 2(- 1 15 a function of x that will satisfy Legendre's 
equation. 

,$` ])2 — ע2)2 — 1[ = )2+ 212 — 1( 


and can be developed by the Binomial Theorem; the coefficient 
of 2" is easily picked out, and proves to be precisely the func- 
tion of x which in Art. 4 we have represented by £,, (+), and 
have called a Surface Zonal Harmonic. 

We have, then, 


)1 -- 222 --2'(-1 =P AHP) .2-- Pl) . ¿EPIA KR... (5) 
if the absolute value of 2 is less than I. 


If x = 1, (5) reduces to 
(1 — 22+ 2'(7-1 = 2%)1(--2)1(.2-7- AD 2+ Pr). 2+...; 
but )1 - 24+ 71=)1(י2‎ - (7 =1 4244284... 
hence Ps (1) = 1. (6) 


Any Surface Zonal Harmonic may be obtained from the 
two of next lower orders by the aid of the formula 


DxP, (1) + 1P,_(x)=0, (7)‏ + א2) -- )3( , 1(15 - א) 


which 15 easily obtained, and 1s convenient when the numerical 
value of x is given. 


Differentiate (5) with respect to 2, and we get 


Say = PH) + 2P (2) 2+ BP) EEE 
whence 
/ Eh —212 4+2 APH 2P r). 243P enn), 


or by (5) 
(1 — 2x2 + PP (4) + 2P(x).24+ 3P) z +.“ 
+ (2 — PA + P). 2 + Pl). + eemo (8) 
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Now (8) is identically true, hence the coefficient of each 
power of 2 must vanish. Picking out the coefficient of 2" and 
writing it equal to zero, we have formula (7) above. 

By the aid of (7) a table of Zonal Harmonics is easily com- 
puted since we have P,(x) = 1, and P(x) = x Such a table 
for x = cos @ is given at the end of this chapter. 


ART. 13. PROBLEMS IN ZONAL HARMONICS. 


In any problem on Potential if V is independent of œ so 
that we can use the form of Laplace’s Equation employed in 
Art. 4, and if the value of V on the axis of X is known, and 


\ m . 
can be expressed as 20," or as > ונק‎ WE Can write out 
‘our required solution as 


V = 3a,r”P,(cosó) or V = > 7 dit add 


--- 4 ) 

for since P,,(1) = I each of these forms reduces to the proper 
value on the axis; and as we have seen in Art. 4 each of them 
satisfies the reduced form of Laplace's Equation. 


As an example, let us suppose a statical charge of M units 
‘of electricity placed on a conductor in the form of a thin circu- 
lar disk, and let it be required to find the value of the Poten- 
tial Function at any point in the “ field of force” due to the 
charge. 

The surface density at a point of the plate at a distance 7 
from its centre is 


M 
c= > 
Aan Va — r 

| TM a 
and all points of the conductor are at potential os See Pierce’s 
Newtonian Potential Function (§ 61). 


The value of the potential function at a point in the axis 
01.66 plate at the distance x from the plate can be obtained 
without difficulty by a simple integration, and proves to be 


(1) 
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The second member of (1) is easily developed into a power 


series. 
A cos-1% e 
24 2% -ך‎ a 
M| 7 x r? x x l 
ב "ל די‎ ca) og . | itz <a (2) 
a @ a a ו‎ 
alr 3r cat ע )1 יה על‎ < 4. (3) 
Hence 


Mix y I 7° 
we 86% 5 | | 
5 q 6 % 0 9 | (4) 


@ 9 o o » T o 
is our required solution if 7 <a and 6 > 57 as 1s 


Mla I 


a Ia 
- — P 6 ו‎ 
"יקי גו‎ (cos 6) + aP P, (cos 0) 


- Pecos קיית‎ 1 ifr >a. (5) 


The series in (4) and (5) are convergent, since they may be 
obtained from the convergent series (2) and (3) by multiplying 
the terms by a set of quantities no one of which exceeds one 
in absolute value. For it will be shown in the next article that 
Pa (cos 0) always lies between I and — 1. 


Prob. 19. Find the value of the Potential Function due to the 
attraction of a material circular ring of small cross-section. 


The value on the axis of the ring can be obtained by a simple 


M 
integration, and is לקו‎ if AZ is the mass and ¢ the radius of the 
c +r 


ring. .At any point in space, 1] 7 > 0 
M Ir 1.3 y 
4 | 0( 7 7 ,)605 0) + שאה‎ =P (cos 0-... , 


and if r > 4 
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| 5 P,(cos 6) — 1 5P, (cos 0) +3 3 SP (cos 1) - 1 


ART. 14. ADDITIONAL FORMS. 


(a) We have seen in Art. 12 that P,,(+r) is the coefficient of 
2" in the development of (1 — 2x2 + 2'(- 1 in a power series. 


Q —212 + 2)-'= [1 — 260% + e7’ + 2’)-3 
= (1 — 2e%)- 1 — ze—*) + 


If we develop (1 — 26"(-+ and (1 — ze-%)-4 by the Bi- 
nomial Theorem their product will give a development for 
(1 — 2xz-+2°)-4. The coefficient of 2" is easily picked out 
and reduced, and we get 


Ps, (COS 0) == 


135 2 605 7/0 + 2 l. 7% 


4.6... Zm T. m — 1) sm 2)0 


I. 3. m(m — 1) 


-++ 2 1.2.(2m—1)(am — 3) cos (m — 4(6 -- .. | (1) 


If + 7 odd the parenthesis in (1) ends with the term con- 
taining cos 0 ; if m is even, with the term containing cos O, but 
in the latter case the term in question will not be multiplied by 
the factor 2, which is common to all the other terms. 


Since all the coefficients in the second member of (1) are 
positive, P„(cos 0) has its maximum value when 0 = O, and its 
value then has already been shown in Art. 12 to be unity. 
Obviously, then, its minimum value cannot be less than — I. 


(0) If we integrate the value of /,,(%) given in (11), Art. 4, 
m times in succession with respect to xx, the result will be 
1.3.5... (27 — Diz? 
(2m)! 
volving lower powers of x than the mth. 


found to differ from — 1)” by terms in- 


a x*— 1)”. (2) 


Hence’ P,(#) = T בר‎ 
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(c) Other forms for F,, (7), which we give without demon- 
stration, are 


(=r וה ה‎ 
fut) = ml Ba" PIPE (3) 
P(x) = AE + YVx* — 1.605 p]"de. (4) 
2100 4 
Emb) = a) — V’ - 8 pyr 2 


(4) and (5) can be verified without difficulty by expanding 
and integrating. 


ART. 15. DEVELOPMENT IN TERMS OF ZONAL HARMONICS. 


Whenever, as in Art. 4, we have the value of the Potential 
Function given on the surface of a sphere, and this value de- 
pends only on the distance from the extremity of a diameter, 
it becomes necessary to develop a function of O into a series 
of the form 


A,P(cos 0) + A,P,(cos 6( + A,P{cos 6)+...; 


or, what amounts to the same thing, to develop a function of 
x into a series of the form 


A,P(2) + A,P(2) + APO +... 


The problem is entirely analogous to that of development 
in sine-series treated at length in Art. 6, and may be solved by 
the same method. 


Assume f(z) = A,P(2) + 4,200 + 4,0(+... (1) 
for -1 > א‎ > I. Multiply (1) by F,(1)dr and integrate from 
— I to I. We get 


/ ANP, אש(‎ = 4 , / Pr DELI AAA. (2) 


-1 
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We shall show in the next article that 


1 
J P A&P Aed = 0, unless 7/ = 2, 
1 


1 i 

2 
and that [P, (1) da = ------ 
J 211 +- 1 


Hence , ome / HAP wdr. (3) 


It is important to notice here, as in Art. 6, that the method 
we have used in obtaining A,, amounts essentially to deter- 
mining A,,, so that the equation 


Hr) = ALP A) APA) HALA) +e. + ALP) 


shall hold good for n + 1 equidistant values of x between — I 
and I, and taking its limiting value as z is indefinitely in- 
creased. 


ART. 16. FORMULAS FOR DEVELOPMENT. 
We have seen in Art. 4 that z = F,, (xr) is a solution of 
L dre's Equati 0 — x E | + n(n I)e =0. (1 
egendre s Equation — x) mk + I)e = O. (1) 


dP „(x 
E ) | 4 m(m-+ Pal) =0, (2) 


a ל‎ 
Hence dí (I — r’) 


= x 
and “| (1-2 - 1 41 +n(z+1)P,(4) = (3) 


Multiply (2) by P,(+) and (3) by P(x), subtract, transpose, 
and integrate. We have 


[n(m + 1( — n(n + 1([ / P,A2)P,(x)de 


[rere א‎ „dP, PAD Jy 
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- J rere 2) ae (a) 


= 4% dP 
= [PRU = 2) ו(‎ — 2) rl 


yq w(t) AP, (x) 
0 / (1 — dx dx . 


+ / a aP ad) ו‎ de 6) 


by integration by parts, 
= O. 


Hence / ע12/)2(1)2(09‎ = O, )6( 


unless 22 = 2. 


If in (4) we integrate from xv to I instead of from — I to I, 
we get an important formula. 


AP, (x AP, 
1 ao -- 52 1 
/ popada aT ד יד יי‎ A 
and as a special case, since P (x) = 1. 
1 (I — AZ) 
/ Pal (09 — 7 אע‎ + 1) | , )8( 


unless m = O. 


1 
To get /'[P,(x)J'dx is not particularly difficult. By (2), 
ו-‎ 


Art. 14, 
1 1 | 
--- "627 — 1)", du — 1)” 
/ ו‎ => ny de da RN 


By successive integrations by parts, noting that 
שש"‎ 


m, and‏ > א contains (1? — 1(* as a factor if‏ ”)1 — ?1( כ 
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that io al = (277)! we get 
1 o 2 | ; ‘ | 
Jete) pde = + ap J (x? — 1)"dx, (10) 


feiron הלש‎ 


1 
. 2 o fi mol 
הב *י-‎ J 1)” (x + עשי"(‎ 


pmt (772 (( 


= (— 1) 4 Sus 1) "dr = (— EN 1 


oi 


1 
. 0 2 
Hence J [Pal] UN ‘Saas (11) 


-1 


1 
“Prob. 20. Show that |] P,(x)de = o if m is even and is not zero 
0 


14-1 
ו‎ .5.7 - 
= 1 1) 2 - בה‎ + : y a 6 ₪ if 2 is odd. 


1 
è \ I 
Prob. 21. Show that / | 2,,)2( ל = שש[‎ Note that 
0 


[ /2,,)2( |” is an even function of x. 


Prob. 22. Show that if V(x) = 0 from x = — 1 tox = O, and‏ ו 
f(x) = 1 from x= 0 10 x= 1,‏ 


Me) = + 5209 — 709 + a) +. 


Prob. 23. Show that .2)0( = 2 B mP, (cos 0) where 


E ו מ‎ 0) sin 6 dô, 
2 + 
0 
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Prob. 24. Show that 


csc 0 = | + (=) P, (cos 0) + ==) P (cos 6( + .. | 
See (1), Art. 14. | 
Prob. 25. Show that 


n! 


27 I 
x* = קוא קה‎ ree) ₪ + 1)P,,(x) + (22 — 3) Pn-2(x) 
on + 1()28 — 1 ke eae 
+ (en — EDI pia) py... 
: > mfa? 
Note that | )"א‎ = . fa ce DY” dy, and use the 
Y 271 ! 0 ax 


method of integration by parts freely. 


Prob. 26. Show that if V is the value of the Potential Function 
at any point in a field of force, not imbedded in attracting or repel- 
ling matter; and if V = f(0) when 7 = a, 


E SAn Pi Kcos 0) 11 r <a 


and Pa SAP „(cos 0) if 7 > a, 
271 + 1 / | 
where Am = [FOP nl cos 0) sin 020, 
` 0 


Prob. 27. Show that if 
V=c¿whenr=a; V=cifr<a, and y= ifr >a, 


ART. 17. FORMULAS IN. ZONAL HARMONICS. + 


The following formulas which we give without demonstra- 
tion may be found useful for reference: © 


BPM) (an —1) Py. (2) (28 5)P, a2) 420-9) MT 


6 _ oe) = (22 + DP) 


AA E AO 


SPHERICAL HARMONICS, AT 


ART. 18. SPHERICAL HARMONICS. 


In problems in Potential where the value of Y is given on the 
surface of a sphere, but is not independent of the angle fp, we 
have to solve Laplace’s Equation in the form (5), Art. 1, and 
by a treatment analogous to that given in Art. 4 it can be 


proved that 


ae Hi H ; . LP 
ETS cf ja” 


where u = cos 0, are particular solutions of (5), Art. 1 


The factors ‘multiplied by 7” in these values are known as: 
Tesseral Harmonics. They are functions of ø and 6, and they 
play nearly the same part in unsymmetrical problems that the: 
Zonal Harmonics play in those independent of q. 

Vu p) = A,P. pe) + ‘> (A„ cos no +B, sin 1p)sin” EL (H+) 


1/1 = 1 


is known as a Surface Spherical Harmonic of the th degree, 
I 
and PESA p) and V = mp ו‎ P) 


satisfy Laplace’s Equation, (5), Art. I. 


The Tesseral and the Zonal Harmonics are special cases of 
the Spherical Harmonic, as is also a form P,,(cos y) known as 
a Laplace’s Coefficient or a Laplacian: y standing for the angle 
between r and the radius vector 7, of some fixed point. 

For the properties and uses of Spherical Harmonics we 
refer the student to more extended treatises, namely, to 
Ferrers Spherical Harmonics, to Heine's Kugelfunctionen, or 
to Byerly’s Fourier's Series and Spherical Harmonics. 


ART. 19.* BESSEL’S FUNCTIONS. PROPERTIES. 


We have seen in Art. 5 that 2 = / (x) where 


xv x’ x 
x)= I — — + -— — ——— +... | 
VA ) 2? . 22.4 2.4 .6' | )1( 


* The student should review Art. 5 before reading this article. 
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is a solution of the equation 

d'z | 1 az 
oe == 2 
e 2) 


and we have called J (x) a Bessel's Function or Cylindrical 


Harmonic of the zero order. 


wt) = x 4 
/,)2( = = ==} רמו‎ ee מ שר‎ WI — a Fo. JO 


is called a Bessel’s Function of the first order, and 
z =J (1) 
is a solution of the equation 


d'z 1 dz I\, 
Geta ae tg =O (4) 


x dx 


which is the result of differentiating (2) with respect to x. 


A table giving values of /,(x) and J (r) will be found at 
the end of this chapter. 

If we write /,(x) for z in equation (2), then multiply 
through by raz and integrate from zero to v, simplifying the 
resulting equation by אוו‎ by parts, we get 


or, since J (1) = = m 
JxI kede = zJ (2). (5) 


If we write hol for 2 in equation (2), then multiply through 


=, 


by x from zero to vx, simplifying by inte-- 


pus parts, we get 


A + Cue] — foem 


o = y. © 
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If we replace x by ux in (2) it becomes 


d'z Idz o 
TS — -ד‎ #z=0 (7) 


xdx 
(See (8), Art. .(ף‎ Hencez = Jux) is a solution of (7). 
If we substitute in turn in (7) (uxx) and J (ux) for z, mul. 
tiply the first equation by x/,(1,1), the second by x/,(M,x), 


subtract the second from the first, simplify by integration by 
parts, and reduce, we get 


/ aJ (Ht Y (MX) ax 


=> lead ADAM (Mea (Ha) (8) 


Hence if 4, and m, are different roots of / (ma) = O, or of 
J ua) =0, or of paf (ua) — ÀJ (ma) = o, 


/ AD Mat) J (ux)dr = 0. (9) 


We give without demonstration the following formulas, 
which are sometimes useful: 


kx) = aJ 00s (x cos p)do. (10) 
J (4) = = f sin’ $ cos (x cos pido. (11) 


They can be confirmed by developing cos (x cos ¿), inte- 
grating, and comparing with (1) and (3). 


ART. 20. APPLICATIONS OF BESSEL’S FUNCTIONS. 


(a) The problem of Art. § is a special case of the following: 
The convex surface and one base of a cylinder of radius «a 
and length 4 are kept at the constant temperature zero, the 
temperature at each point of the other base is a given function 
of the distante of the point from the center of the base; re- 
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quired the temperature of any point of the cylinder after the 
permanent temperatures.have been established. 


Here we have to solve Laplace’s Equation in the form 


O` I Ou | Ou _ 
art par on 9 1) 


(see Art. 5), subject to the conditions 


u =o when z=0, 

4=0 * FH a, 

/ = Ar)“ 2=0. 
Starting with the particular solution of (1), 


u = sinh (uz) / (ur), (2) 
and proceeding as in Art. 5, we get, if U,, Mz, Mas 2. are roots 
of J (Ma) = O, (3) 


and f(r) = 20, Jdr) + AS lu) + AS Mr) +--+, (4) 


sinh (4,2) sinh (1, 2) 
ETICO יי‎ ‘sinh (u,b) וו‎ (5) 
(9) If instead of keeping the convex surface of the cylinder 
at temperature zero we surround it by a jacket impervious to 
heat the equation of condition, א‎ = O when 7 = a, will be re- 


placed by a= o when r = a, or if u = sinh (yz)/,(uz) by 


aS (Ur) =0 whenv=a, 
ar 
that is, by — p/ fa) =0 
or J (Ha) = 0. (6) 


If now in (4) and (5) M,» AM,» Ms, - - - are roots of (6), (5) will 
be the solution of our new problem. 


(c) If instead of keeping the convex surface of the cylinder 
at the temperature zero we allow it to cool in air which is at 
the temperature zero, the condition u =0 when z =a will be 


replaced by = + hu = 0 when r = a, being the coefficient 


of surface conductivity. 
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If u = sinh (uz)/,(u7) this condition becomes 
— MJ (ur) + A/ (Mr) = 0 when ? =a, 
or paj (ua) — aj (ua) = 0 (7) 
If now in (4) and (5) 4,, Ma, My, - - - are roots of (7), (5) will 


be the solution of our present problem. 
It can be shown that 


J (4) = 9, (8) 
J (4) = 9, (9) 
and | a] (1) - (ע),/4‎ = 0 )10( 


have each an infinite number of real positive roots.* The 
earlier roots of these equations can be obtained without serious 
difficulty from the table for /,(7) and / (r) at the end of this 
chapter. 


ART. 21. DEVELOPMENT IN TERMS OF BESSEL'S FUNCTIONS. 


We shall now obtain the developments called for in the last 
article. 


Let Ar) =A Lar) + AL (147) + ASM) Fee מ‎ 
H, Ma» Hz, etc., being roots of / (ua) = o, or of / (ua) = O, or 
of paj (ua) — AJ (ua) = 0. 


To determine any coefficient 4, multiply (1) by ~/7,(4u7)dr 
and integrate from zero to a. The first member will become 


[fT und 


0 
Every term of the second member will vanish by (9), Art. 
10, except the term 


A, fr [J (mr) dr. 


JT = zz fas inter =F (ua +a?) 
by (6), Art. 19. | 


* See Riemann’s פיפ‎ 86 Differentialgleichungen, S 97» 
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2 a 
— | f(r iv) ar. 
וקו‎ aro A Jr Arde. (2) 


The development (1) holds good from r = 0 to r = a (see 
Arts. 6 and 15). 


If M, My» Ms, etc., are roots of / (ua) = 0, (2) reduces to 


Hence A, = 


2 a 
A, = ה‎ rfr) (Mar ar | )3( 
If p,, ,מ‎ etc., are roots of / (ua) = O, (2) reduces to 
A, = A [fry ) 0707 (4) 
T Gaf) SVMs 


If Mis My, My, etc., are roots of paf (ua) — ÀJ (ua) = O, 
(2) reduces to ' 


POE? __ 
=p) TO 9 


For the important case where /(7) = 1 


Mpa 


/ AS, quoyar= fr, ear a fad (1)dr = Ta) (6) 


by (5), Art. 19; and (3) reduces to 


2 
A, == — r; 
paa] (usa) 0( 
(4) reduces to 
A, = O, (8) 
-except for א‎ = I, when u, = O, and we have 
A, = A. (9) 
24 


(5) reduces to (10) 


4 = + uay ma) 
Prob. 28. A cylinder of radius one meter and altitude one metër 
has its upper surface kept at the temperature ,"ססז‎ and its base and 
convex surface at the temperature 15°, until the יי‎ temper- 
atures are established. Find the temperature at points on the axis 
25. 50, and 75 centimeters from the base, and also at a point 25 
.centimeters from the base and 50 centimeters from the axis. 
j Ans. 29.6; 47°.63 7172 ; 8 
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Prob. 29. An iron cylinder one meter long and 20 centimeters 
an diameter has its convex surface covered with a so-called non-con- 
ducting cement one centimeter thick. One end and the convex 
surface of the cylinder thus coated are kept at the temperature zero, 
-the other end at the temperature of roo degrees. Given that the con- 
ductivity of iron 1s 0.185 and of cement 0.000162 in C. G. S. units. 

Find to the nearest tenth of a degree the temperature of the mid- 
dle point of the axis, and of the points of the axis 20 centimeters 
from each end after the temperatures have ceased to change. 

Find also the temperature of a point on the surface midway be- 
tween the ends, and of points of the surface 20 centimeters from 
each end. Find the temperatures of the three points of the axis, 
supposing the coating a perfect non-conductor, and again, suppos- 
ing the coating absent. Neglect the curvature of the coating. Ans. 


15.453 40.85; 72.8; 15.33 40.73 72-53 0-03; 0.0; 1.3. 
Prob. 30. If the temperature at any point in an infinitely long 

cylinder of radius ¢ 15 initially a function of the distance of the 

point from the axis, the temperature at any time must satisfy the 


Ou, I Ou 


equation 0 a | (see Art. 1), since it is clearly in- 


7 or’ * yor 
dependent of z and Q. 
Show that 


u = Ae] (ar) + Aye vet] (uy) 
+ Aj UU +... 


where, if the surface of the cylinder is kept at the temperature 
zero, U,, M,, U,,... are roots of J, (uc) = o and A; is the value 
‘given in (3) with ¢ written in place of a; if the surface of the cylin- 
der is adiabatic M,» M,, M,, » -- are roots of /, (Mc) = o and 4; is ob- 
tained from (4); and if heat escapes at the surface into air at the tem- 
perature zero 4,» Hy, Haee. are roots of Mc) (Mc) — AJ (pc) = o, 
and A; is obtained from (5). 


Prob. 31. If the cylinder described in problem 29 is very long 
and is initially at the temperature 100” throughout, and the con- 
vex surface is kept at the temperature o”, find the temperature of a 
point 5 centimeters from the axis 15 minutes after cooling has begun ; 
first when the cylinder is coated, and second, when the coating 5 
absent. Ans. 07%.2 ; 1. 

Prob. 32. A circular drumhead of radius æ is initially slightly 
‘distorted into a given form which 15 a surface of revolution about 
the axis of the drum, and is then allowed to vibrate, and 2 is the 
ordinate of any point of the membrane at any time. Assuming that 
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2 


₪ 
2 must satisfy the equation = z 


art 
- OZ | 
ditions 2 = o whens = a, 5, = 0 when / = o, and z = f(r) when: 
t = o, show that z= 4A, J (ur) cos uct + A, J (Har) cos uet +.. 
where 4,, M,, M,, +. are roots of /, (Ma) = o and 4; has the value 
given in (3). 

Prob. 33. Show that if a drumhead be initially distorted as in 
problem 32 it will not in general give a musical note ; that it may be 
initially distorted so as to give a musical note ; that in this case the 
vibration will be a steady vibration ; that the periods of the various 
musical notes that can be given are proportional to the roots of 
J (x) = 0, and that the possible nodal lines for such vibrations 
are concentric circles whose radii are proportional to the roots of 


Jo (x) = ₪ 


ART. 22. PROBLEMS IN BESSEL’S FUNCTIONS. 


ht , 4 subject to the con--‏ א 


If ina problem on the stationary temperatures of a cylinder 
u = 0 when 2 = 0, א‎ = O when 2 = J, and u = f (2) when 7 = a, 
the problem is easily solved. If in (2), Art. 20, and in the cor- 
responding solution 2 = cosh (vz)/,(u7) we replace u by uz, we 
can readily obtain z= sin (uz)/,(ur2) and 2 = cos (12)/,(u71) 
as particular solutions of - Art. 20; and 


| x 
A(z) = ג‎ +> + + soe ווה‎ (1} 
and is real. 
f(z) = = Ax ₪ יש"‎ 
2 l RTZ | 
where Ay =z JF sin גש‎ dz 12) 
by Art. Q. 
Rar? 
Luo. ו‎ bag 0 0 | 
Hence u = = A, Sin "+" 1-2 13) 
Ja ₪ 


is the required solution. 
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A table giving the values of /,(x2) will be found at the end 
of this chapter. 

Prob. 34. A cylinder two feet long and two feet in diameter has 
its bases kept at the temperature zero and its convex surface at 
100 degrees Centigrade until the internal temperatures have ceased 
to change. Find the temperature of a point on the axis half way 
between the bases, and of a point six inches from the axis, half way 
between the bases. Ans. 72. 1; 80°.1. 


ART. 23. BESSEL'S FUNCTIONS OF HIGHER ORDER. 


If we are dealing with Laplace's Equation in Cylindrical 
‘Coordinates and the problem is not symmetrical about an 
axis, functions of the form 


A A A 
oe 2"] (z+ 1) 2% -- 1)" 2.260 Hia) 
play very much the same part as that played by /,(x) in the 
preceding articles. They are known as 2655613 Functions of 
the mth order. In problems concerning hollow cylinders much 
‘more complicated functions enter, known as Bessel's Functions 
of the second kind. 

For a very brief discussion of these functions the reader is 
referred to Byerly’s Fourier's Series and Spherical Harmonics ; 
for a much more complete treatment to Gray and Matthews’ 
admirable treatise on Bessel’s Functions. 


ART. 24. LAME’S FUNCTIONS. 


Complicated problems in Potential and in allied subjects are 
usually handled by the aid of various forms of curvilinear co- 
órdinates, and each form has its appropriate Harmonic Func- 
> tions, which are usually extremely complicated. For instance, 
Lamés Functions or Ellipsoidal Harmonics are used when 
solutions of Laplace’s Equatiomin Ellipsoidal coordinates are 
required; Toroidal Harmonics when solutions of Laplace’s 
Equation in Toroidal coordinates are needed. 


For a brief introduction to the theory of these functions 
566 Byerly's Fouriér’s Series and Spherical Harmonics. 
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TABLE I. SURFACE ZONAL HARMONICS. 

P, (cos 6(| Py cos 8) | Py (cos 0) P, (cos 6) | P, (cos 8) P, (cos 0) P, (cos 0) 
7071 | .2500| —.1768 | —.4062 | —.3757 | —.1485 1270 
6947 | .2238 | —.2040 | --.4158 | —.3568 | —.1079 1666 
6820 | .1977 | --.9900 | —.4252 | —.3350 | —.0645 2054 
6691 | .1716 | —.2547 | —.4270 | —.3105 | —.0251 2349 
6561 | .1456 | —.2781 | —.4286 | —.2836 0161 2627 
.6428 | .1198 | —.3002 | —.4275 | —.2545 0563 2854 
6293 | .0941 | —.3209 | —.4239 | —.2235 0954 3031 
.6157 | .0686 | —.3401 | —.4178 | —.1910 1326 8153 
.6018 | .0433 | —.3578 | —.4093 | —.1571 1677 3221 
.5878 | .0182 | --.8740 | —.3984 | —.1223 2002 3234 
.5736 | —.0065 | --.8886 | —.3852 | —.0868 2297 8191 
.5595 | —.0310 | —.4016 | —.3698 0510 2559 8095. 
5446 | —.0551 | -.4191 | —.3524 | .0 2787 2949 
5299 | --.07-8 | —.4229 | —.3331 0206 2976 752 
.5150 | —.1021 | —.4310 | —.3119 0557 8125 511 
5000 | —.1250 | —.4375 | —.2891 0898 8232 2931 
4848 | --.1474 | —.4423 | —.2647 1229 8298 1916 
4695 | —.1694 ' —.4455 | —.2390 1545 3321 1 
4510 | —.1908 —.4471 | —.2121 1844 3302 1208 
.4884 | —.2117 | —.4470 | —.1841 2123 3240 0818 

| 
‘4226 | —.2321 | —.4452 | —.1552 2381 3188 0422 
.4067 | —.2518 | —.4419 | —.1256 2615 .2996 0021 
3907 | —.2710 | —.4370 | —.0955 2824 2819 | .3 
.3746 | —.2896 | —.4305 | —.0650 3005 2605 | —.0763 
.3584 | —.3074 | —.4225 | —.0344 3158 2361 | —.1135 
.9490 | —.3245 | —.4130 | —.0088 8281 2089 | —.1485 
.3256 | —.3410 — 14021 | 0267 .3373 1786 | —.1811 
.3090 | —.3568 | —.3898 0568 8434 1472 | .9 
2924 | —.3718 | —.3761 4 .8 1144 | —.2347 
.2756 | —.3860 | —.3611 1153 .3461 0765 | —.2559 
2588 | —.3995 | —.3449 1434 8427 0481 | —.2730 
2419 | —.4112 | —.3275 1705 3362 0076 | —.2848 
2250 | —.4241 | —.3090 1964 8267 | —.0284 | .9 
2079 | —.4352 | —.2894| .2211 8143 | —.0644 | —.2943 
.1908 | —.4454 | --.9688 | 9 .2990 | —.0989 | .8 
.1736 | --.4548 | —.2474|  .2659 2810 | —.1321 | .% 
.1564 | —.4633 | - .1 9 2606 | —.1635 | —.2709 
.1392 | —.4709 | —.2020 8040 2378 | —.1926 | —.2536 
1219 | —.4777 | —.1783| 8 .2129 | —.2193 | —.2321 
1045 | —.4836 | —.1539 8345 1861 | —.2431 | —.2067 
0872 | —.4886 | —.1291 |  .3468 1577 | —.2638 | —.1779 
.0698 | —.4927 | 8 .3569 1278 | —.2811 | —.1460 
0523 | —.4959 | —.0781 3648 0969 | —.2947 | —.1117 
.0949 | —.4982 | —.0522 3704 0651 | —.3045 | —.0735 
.0175 | —.4995 | —.0262 3739 .0397 | —.3105 | —.0381 
.0000 | —.5000 | .0 .3750 .0000 | —.3125 0000 
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TABLE II. BESSEL’s FUNCTIONS. 


x J (2) NET x J o(X) J (x) x | J (Xx) Jia) 

0.0 | 1.0000| 0.0000 | 5.0 | —.1776 | —.3276 | 10.0 | —.2459 | 8 
0.1 9975} .0499 | 5.1 | —.1443 | —.3871 | 10.1 ' —.2490!  .0184 
0.2 .9900 | .0995 | 5.2 | —.1103 | —.3432 | 10.2 | —.2496 | .0066 
0.3 9776! .1483 | 5.3 | --.0788 | —.3460 | 10.3 | —.2477 | —.0813 
0.4 .9604 | | .1960 | 5.4 | —.0412 | —.3453 || 10.4 | —.2434 | —.0555 
0.5 .9385 | .2423 | 5.5 | --.0068 | —.3414 |i 10.5 | —.2366 | —.0789 
0.6 9120 | .2867 | 5.6 | .0270| —.3343| 10.6 | —.2276 | —.1012 
0.7 .8812 | .3290 | 5.7 | .0599| —.3241 | 10.7 | --. 1664: .4 
0.8 .8463 .3688 9.8 -0917) —.3110 | 10.8 | —.2032 | —,1422 
0.9 .8075 | .4060 | 5.9 | .1220| —.2951 |; 10.9 | —.1881 | —.1604 
1.0 | .7652} .4401 | "6.0 | .1506| —.2767| 11.0 | —.1712 | —.1768 
1.1 7196 | 14709 | 61 | /1773|—.2559| 11.1 | —.1528 | —.1918 
1.2 .6711| .4983 | 6.2 | .2017| —.2329 | 11.2 | —.1830 | —.2039 
1.3 6201 | .5220 | 6.3 .2238 | —.2081 | 11.8 | —.1121 | —.2143 
1.4 .8669 | .5419 | 6.4 | .2433| --.1816 | 11.4 | —.0902 | —.2225 
1.5 .5118 | .5579 | 6.5 | .2601| —.1538| 11.5 | —.0677 | —.2284 
1.6 4554| 5699 | 6.6 | .2740|—.1250| 11.6 | —.0446 | —.2320 
14 .3980 | | .5778 | 6.7 | .2851 | —.0953 || 11.7 | —.0213 | .8 
18 | .8400| .5815 | 6.8 | 2931 | —.0652 | 11.8 | 0 8 
1.9 .9818 | .5812 | 69 | .2981 | —.0349 | 11.9 | .0250 | —.2290 
2.0 2239 | .5767 | 7.0 | .3001 | --.0047 | 120 |  .0477 | —.£284 
1 .1666 | .5683 | 7.1 | .2991| .0252, 12.1 | .0697 —.2157 
2.2 .1104| .5560; 7.2 | ' .2951 | .0543: 122 | 0908! —.2160 
2.3 .0555 | .5899 | 7.3 | .2882| .0826! 199 | .1108| .8 
4 .0025 | .5202 | 7.4 | .2586| .1096| 4 1296 —.1807 
95 | —.0484| .4971 | 7.5 | .2668| .1352| 19.5 | .1469| —.1655 
2.6 | —.0968| .4708 7.6 | .2516| .1592| 12.6 1620 —.1487 
27 | -.1424| .4416 | 7.7 | .2346| .1813| 127 | .1766 | —.1307 
98 | —.1850| .4097 | 7.8 | .2154| .2014¡ 128 | .1887 | —.1114 
2.9 | -.9948 | .3754 | 7.9 | .1944| .2192| 129] .1988| .2 
8.0 | -.9601| .8891 | 80 | .1717| .2346| 19.0 | .9069| .8 
31 | -.9991| .3009 | 8.1 | .1475| .2476' 181 | .21299| — .9 
3.2 | -.3909 | .2613 | 89 | .1222| .2580| 18.9 | 9167 | .1 
33 | -.9443 | .2207| 8.3 | .0960| .2657 | 18.8 | .2188| —.0052 
3.4 | —.3643/ .1792| 84 | .0692| .2708 || 13.4 ed .0166 
35 | —.3801| .1374|| 8.5 | .0419| .2781 || 135 | .2150| .0880 
3.6 | —.3918| .0955|| 8.6 | .0146| .2728¡| 18.6 | .2101| 0 
3.7 | —.3992| .0538| 8.7 | -.0195 | .2697 || 13.7 | .2082] -.0791 
3.8 | —.4026 | .0128 || 88 | -.0892 | .2641 | 188 | .1948]  .0984 
39 | —.4018| —.0272 | 8.9 | —.0653 2559 | 13.9 | .1886| .6 
40 | -.8979 | —.0660 || 9.0 | -.0908 | .2458'| 140 | .1711| .4 
4.1 | —.3887| —.1033 || 9.1 | —.1142] .23894 (| 111 | .1570| .8 
4.2 | —.3766| —.1386 | 92 | -.1967 | .2174 i| 149 | .1414| .86 
4.3 | —.3610| -.1719 || 9.3 | —.1577| .2004 | 14.8 | .1245| .1747 
4.4 | —.3423| —.2028 || 9.4 | -.1768 | .1816 || 14.4] .1065| 0 
4.5 | —.3205| —.2311 || 9.5 | —.1939| .1613 || 145 | .0875| .1984 
4.6 | —.2961 | —.2566 || 9.6 | - .9090 | .1895! 14.6 | .0679| .1999 
4.7 | -.9698| —.2791 | 9.7 | —.2218| .1166 || 147 | .0476|  .2043 
4.8 | —.2404| --.9985 | 9.8 | —.2323| .0928 || 14.8 | .0271| 6 
4.9 | —.2097 | —.3147 || 9.9 | -.2403 | .0684 ¡| 14.9 | .0064| 9 
5.0 | —.1776 | —.3276 || 10.0 | -.9459 | .0485 || 15.0 | --.0149| 1 


TABLES. 03. 


TABLE III.—Roots OF BEssEL’s FUNCTIONS. 


| 
n x, for Jo(X,) = 0 |x, for J\(x,) = 0 n x, for Jo(x,) = 0 |x, for J¡(x,) = 0 


1 2.4048 3 8317 6 18.0711 19.6159 
2 9.9201 . 6 Y 21.2116 22. (601 
5 8.6531 10.1735 8 24.3025 20 9087 
4 11.7915 13.3237 9 21.4935 29.0468 
0 14.9309 : 16.4706 10 80.6346 ox. 1897 


TABLE IV.—VALUES OF Jo(zi). 


0.0 1.0000 2.0 2.2196 4.0 11.3019 
0.1 1.0025 2.1 2.4463 4.1 12.3236 
0.2 1.0100 2.2 2.6291 4.2 13.4425 
0.3 1.0226 2.0 2.8296 4.3 14.6680 
Q 4 1.0404 2.4 9.0493 4.4 16.0104 
0.3. 1.0635 2.0 3.2898 4.5 17.4812 
0.6 1.0920 2.6 3.5993 4.6 19.0926 
0 7 1.1263 2.7 9.8417 4.7 20.8585 
0.8 1.1665 2.8 4 3 4.8 22.7937 
0.9 1.2130 2.9 4.5027 4.9 24.9148 
1.0 1.2661 3.0 4.8808 5.0 27.2399 
1.1 1.3262 3.1 9.2945 5.1 29.17889 
1.2 1.3987 3.2 9 7472 5.2 99 . 0836 
1.3 1.4963 3.9 6.2426 9.8 30.6481 
14 1.5584 3.4 6.7848 9.4 39.0088 
1.5 1.6467 3.0 7.3782 5.5 42.6946 
1.6 1.7500 3.6 8.0277 5.6 46 .6 
1.7 1.8640 25.7 8.7386 D.T 1.1725 
1.8 1.9896 3.8 9.5169 5.8 06.0381 
1.9 2.1277 3.9 10.3690 9.9 61.3766 
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zonal, 12, 15-21, 40, 50, 60-61. 
Heat v. Conduction of heat, 7 
Historical introduction, 7. 


Introduction, historical and descriptive, 
7, 8, 9- 


Lamé, 7. 
Lamé’s functions, 12, 59. 
Laplace, y. 


66 INDEX. 


Laplace’s coefficients, 12, 1. 
Laplace’s equation, 17, 41, 43, 5I. 
in cylindrical coordinates, 10, 21. 
in spherical coordinates, 9, 12. 
Laplacian, 51. 
Legendre, 7. 
Legendre’s coefficients, 19. 
Legendre’s equation, 17, 40, 41, 47. 
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Stationary temperatures: 
problems, 21, 25, 56, 57, 59. 
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